Effect of dielectric responses on localization in ID random periodic-on-average systems 
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Dielectric response effects on wave localization in random periodic-on-average layered systems 
(POAS) are studied. Based on Monte Carlo simulations and products of Random Matrices, statis- 
tics of the Lyapunov exponent are determined efficiently for very long systems. A novel oscillatory 
behavior for Lyapunov exponent is found and explained for mildly strong scattering conditions. We 
also show the emergence of strongly localized states in metallic layered systems with intermediate 
disorder for frequencies above the plasma frequency u) p of metals, as is not shown in dielectrics. 
Furthermore, the violation of universal single parameter scaling behaviors in different regimes of 
multiple scattering is discussed. 
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The realization that Anderson localization in electronic 
systemsEl is due to wave interference between scattered 
waves from random scatters has stimulated vivid re- 
search in search of wave localization in condensed mat- 
ter physicseTEl. Further progress-have been made to the 
aspects of photon localizationEHiJ because photons may 
open a new realm of optical and microwave phenomena, 
and provided an analogous insight into Anderson local- 
ization transition undisturbed by Coulomb interaction. 

It had been shown individual scatter's response to the 
wave fields influence wave localization properties remark- 
ably in entirely random systems , such as the_strong local- 
ization for acoustic waves in bubbly liquids!! The same 
issues should exist in POAS for the intricate interplay 
between order and disorder. One-dimensional systems 
are particular interesting as they provide insights to the 
problems of wave localization in general and are suit- 
able for testing various ideas. Two qualitatively different 
regimes of localization are exhibitedOEJ. For band gap 
states, single parameter scaling (SPS) with universal be- 
haviors is observed; however, the scaling is restored only 
when the randomness of defects exceeds a certain thresh- 
old for the situation of weak dielectric mismatch between 
constituent layersta. 

The scope of the Communication is twofold. The main 
interest is to understand multiple scattering effects on 
wave localization and the SPS in ID systems. Particular 
attentions are paid for mildly scattering conditions; that 
is, the dielectric mismatch between layers is not too large. 
The secondary interest lies in preliminary exploration 
of dispersive or absorptive media on wave localization. 
Here, we consider EM wave localization in ID random su- 
perlattices made of two alternating layers, the so-called 
Kronig-Penny model. As usual, the wave transmission 
can be .tackled in exact manner by the transfer-matrix 
methodEJ. Nevertheless, the approach encounter the dif- 
ficulty that transmission coefficient in that formulation 
falls bellow the computer round-off accuracy for very long 



systems 



For the practical purpose, we improve the 



Herbert-Jones- Thouless formulaO widely known in the 
study of disordered electronic quantum systems to study 
electromagnetic wave localization, with electron energies 
being replaced by wave frequencies; consequently, uni- 
versal behaviors for very long structures can be studied 
quantitatively. Numerically, the state of wave fields is 
formulated by a vector v n = (E n ,E n y with the elec- 
tric field E n , its spatial derivative E n , the wave number 
k = u)/c in vacuum, and the superscript t represents the 
matrix transposition. The relation between the (n+l)th 
layer and the nth layer is specified by the transfer matrix 
M n : 



M,, 



cos(k n d n ) £-sin(fc n d„) 
-k n sm(k n d n ) cos(k n d n ) 



(1) 



where k n — k\J e n (u>), and the n-th layer width d n . 
In dielectrics, the dielectric constant e n is real; how- 
ever, in metals, the dielectric constant e n is complex and 
frequency-dependent. Lyapunov exponent (LE) and its 
variance (VAR), Var(A) = (A 2 ) — (A) 2 , are the transmis- 
sion quantities considered here. The VAR represents the 
fluctuations for the LE in defect configurations. The LE 
can be computed as followstj: 



2N 



A = Re lim 

N-«x: 



In Tr J| M„(Conj(e„)) 



(2) 



n=l 



where 2N is the total number of alternating layers, L is 
the average total length of considered systems, Tr is the 
trace operator, Conj means taking complex conjugate, 
and the symbol (...) means the ensemble average over 
the uniform distribution of random configurations. The 
complex conjugate of the dispersive dielectric function_e„ 
is taken for the reciprocity of the absorptive systemsEj. 

Consider ID superlattices made of two alternating lay- 
ers A and B with the dielectric response functions ca and 
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€b respectively. Random configurations are introduced 
by varying the layer B width ds uniformly in the interval 
({de)(l — &), (^b)(1 + o)) , where a G (0, 1) represents 
the randomness while the width of A layer <1a is kept 
fixed. A series of Random Matrices M„ are generated 
to obtain the LE and VAR by the Monte Carlo simula- 
tion in Eq.(2). In practice, the accurate LE and VAR 
are attained by letting the total length of layered sys- 
tems to be several thousand times larger than the decay- 
ing length-scale A -1 (the inverse of Lyapunov exponent), 
representing exactly the localization length in absorption 
free cases. Sufficient times of ensemble averages are thus 
determined by the evaluation of numerical stability with 
the length of the system. 
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FIG. 1. Model 1: (al) Band structures where K is Bloch 
wave number in periodic layers and k = u/c is the wave num- 
ber in free space, (bl) Lyapunov exponent vs. randomness 
a. (cl) Variance vs. randomness a. Model 2 : (a2) Band 
Structures. (b2) Lyapunov exponent vs. randomness a. (c2) 
Variance vs. a. 

To begin with, multiple scattering effects on wave lo- 
calization in dielectrics are studied. Two models are 
demonstrated. Model 1 is a weakly scattering model 
with the dielectric contrast es/^A — L21 while Model 
2 is the mildly strong scattering model with the contrast 
e^/eA- j= 6.25. The band structures of original periodic 
layersEJ for Model 1 and 2 are shown in Figs. 1 (al) 
and (a2). The LE and VAR in Model 1 and 2 are com- 



pared in Figs. l(bl - cl) and (b2 - c2). The parameters, 
6a = l,d,A = ((1b) = d = 1, are taken in computations. 
The symbols V,'o', and 'o' denote states in bands, band 
edges, and gaps respectively hereafter. .— . 

In Model 1, previous results are reproducedli-l In Fig. 
1 (bl), the LE decreases monotonically with randomness 
in the first gap around k = 1.50 since the established 
Bragg wave coherence has been gradually destcoyed, the 
so-called enhanced transmission due to disorder!!!]. In the 
allowed band k = 1.43, the LE initially increases with the 
randomness and then decreases for large randomness. In 
the band edge k = 1.454, the LE decreases first, then 
increase slightly, and finally decrease monotonically with 
the randomness a. The behavior is due to the interplay 
between order and disorder in the systems. Moreover, 
the curves for these states merge for large disorder- the 
manifestation of universal single parameter scalingii In 
Fig. 1 (cl), the plot of VAR vs. randomness a is demon- 
strated. The VAR at these frequencies shows similar be- 
haviors. When the randomness of defects is large in the 
ensembles, the effect of self-averaging is effective because 
random fluctuations for these states are almost cancelled 
leading to the reduction in VAR. For nearly periodic con- 
figurations, the wave localization length in a gap regime 
is much smaller than the mean distance between defects 
so that the defects are isolated for incident waves. The 
increase of randomness can only cause minor growth of 
fluctuations due to the increase of the localization length. 
The competition between the two mechanisms results in 
the observed maximal VAR. 




FIG. 2. The average transmission (lnT) = —2\N vs. fre- 
quency k with different randomness a. 

In Model 2, mildly strong scattering case, novel oscil- 
latory behaviors of the LE and VAR in the 3rd gap are 
exhibited above the randomness a = 0.3 in Figs. 1 (b2) 
and (c2), which is very different from the behaviors of 
Model 1 in Figs. 1 (bl) and (cl). To understand the dis- 
order effects, we plot the average transmission spectrum 
(lnT) covering the lowest three pass bands and gaps of 
Model 2 for the randomness a from 0.1 to 1, as is shown 
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in Fig. 2. First, we observe well denned band struc- 
tures, transmission enhancement in gaps, and inhibition 
in pass bands appears as usual when the minor random- 
ness ranging from a = 0.1 to a = 0.3 is presented. With 
large randomness above a — 0.3, we observe the enor- 
mous overlap between original band structures emerges 
from the high frequency states from k — 2 to 3 and then 
extends to the states from k = 1 to 3. In addition, it is 
noticed that the oscillatory behavior only occurs for these 
states. Take the 2nd gap for example, the transmission 
initially enhance till the randomness a — 0.6, and then 
weaken for larger randomness. For the 3rd gap, the oscil- 
lation in the transmission corresponds to the oscillation 
in LE in Fig. 1 (b2). For the first gap and pass band 
where the overlap with other states is negligible, the be- 
havior does not exist. Therefore, the novel oscillation in 
LE could be explained as a result of the overlap between 
band structures in these states for large randomness. 
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FIG. 3. Model 3: (a3) Lyapunov exponent vs. randomness 
a. (b3) Variance vs. randomness a . Model 4: (a4) Lyapunov 
exponent vs. randomness a at the frequency k = 1.630 in 
Model 3. The circle o, triangle A, and square □ represent 
the cases with electron damping rates 7 = 0, 0.02, and 0.04, 
respectively. (b4) Variance vs. randomness a. 

Turn to the secondary interest. As a preliminary in- 
spection on dispersive media, the statistics of LE and 
VAR for metallic slabs are demonstrated in Fig. 3. The 
layers B are composed of metals with the dielectric re- 
sponse function cb '■ 



e B (uj) = 1 



lu(lo + ij) 



(3) 



where uj p is the plasma frequency, ui the EM wave fre- 
quency, and 7 the electron damping rate which leads 
to residual wave energy dissipation. The fixed param- 
eters for layer A are the same as previous Model 1 
and 2. The only difference is the introduction of the 
plasma frequency Lo p d/c = 1 ( c : speed of light in vac- 



uum) of metallic slabs, and the reduced wave number 
k = utd/c = Lu/ujp. 
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FIG. 4. Scaling parameter r vs. randomness for Model 1, 
Model 2, Model 3, and Model 4, respectively. The horizontal 
solid line denotes the universal value r = 1. 

We focus on frequencies above the plasma frequency 
Up since waves can only propagate above the plasma fre- 
quency. It is noted that the dielectric contrast e b I (-a hr 
this case is real in the absence of residual absorption. The 
residual absorption can be turned off by taking electron 
damping rate 7 as zero. The frequencies in later compu- 
tations are taken from the first gap above ui p . In Fig. 3 
(a3), a peculiar feature not shown in dielectrics occurs. 
For the gap state k = 1.634, the contrast between the LE 
at the randomness a ~ and 0.5 is much smaller than 
Model 1 (See Fig. 1 (bl)). Moreover, the LE is even 
larger with the intermediate randomness a = 0.5 than 
the LE with the minor randomness a ~ at the edge 
state k = 1.630. The explanations are as follows. The 
dispersive response of metals makes the dielectric con- 
trast between layers lower for higher frequencies, which 
leads to the reduction of extra transmission difference 
between the pass bands and gaps. So, the feature is at- 
tained by diminishing the difference of these states in LE 
in its dielectric counterpart at minor randomness a ~ 
in qualitative agreement with the behavior in Fig. 3 (a3). 
The effect is remarkable with nearly periodic layers be- 
cause individual layer's responses can be accumulated 
by the Bragg wave coherence. This suggests versatile 
ways could be available in the fabrication of the EM 
wave devices operating near band edges in POAS with 
dispersionsEI In Fig. 3 (b3), the VAR agrees qualita- 
tively with Model 1 although the shape of the VAR is 
not so sharp as Model 1. In addition, the VAR does not 
merge for these frequencies at large degrees of disorder, 
the implication of the SPS violation to be discussed later. 
One may wonder the peculiar feature may be smeared out 
with residual absorption. To confirm the point, the LE 
vs. randomness figure is plotted in Fig. 3 (a4) as com- 
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parison. We see the feature can be strengthened further 
because the LE increases noticeably with the damping 
rate 7 for large randomness. The VAR decreases with 
the dampingp^ate in Fig. 3 (b4), as is shown similarly in 
the referenced . 

Single parameter scaling (SPS) is studied next to un- 
derstand its dependence on constituent scatter proper- 
ties, as is shown in Fig. 4. The SPS indicates that the 
LE is not an independent naxameter and is related to its 
variance in a universal waymLj: 



Var(A) 
A 



L = l, 



(4) 



where L is the length of the system. The horizontal solid 
line represents the universal value r = 1. In Fig. 4, 
for the weakly scattering case, the SPS is obeyed for all 
band gap states in Model 1 . This is reasonable since weak 
gaps make the macroscopic wave coherence_easily broken 
by disorders, as is shown by Deych et alliil However, 
the SPS is not rigorously met for large randomness with 
slight fluctuations in r in Model 2 and Model 3. The phe- 
nomenon is universal in the two models, one is dielectric 
and the other is dispersive, even though their detailed LE 
and VAR disagree. In retrospection, in Model 3, the dis- 
persive dielectric response in metals generates an intrinsi- 
cally larger dielectric mismatch than Model 1. Therefore, 
the degree of deviations near the universal value r = 1 
can be recognized as the signature of the gradual domi- 
nance of multiple scattering effects. In Model 4, we show 
that the parameter r decreases with electron damping 
rate as well as randomness. Moreover, the residual ab- 
sorption makes the value of the parameter r smaller than 
one. From above discussions, the SPS depends much on 
varied dielectric responses of constituents. In weak scat- 
tering cases, the SPS is met for allowed bands, gaps, and 
edges. In mildly scattering cases, the parameter shows 
deviated fluctuations near the universal value r = 1. 
However, in residual absorption cases, the parameter r 
falls below the universal value. 

In addition, for strongly scattering conditions, the SPS 
is violated for two possibilities in our investigations. One 
is due to narrow band structures created in underlying 
periodic layers resulting in tremendous oveclap between 
these bands in the presence of randomnessO. The other 
possibility checked by us, however, is the generation of 
wide bands by exchanging the sequence of the same con- 
stituent layers such that the overlap does not emerge and 
gaps are robust enough to resist the addition of random- 
ness. Therefore, to sum up, the band structures of un- 
derlying periodic layers in various scattering regimes de- 
termines the distinctive wave localization properties in 
POAS. 

In conclusion, we have studied effects of dielectric re- 
sponse on wave localization and single parameter scaling 
behaviors in periodic-on-average systems. Novel oscilla- 
tory wave localization behaviors in the presence of mildly 
strong multiple scattering effects is found and explained 
numerically. For dispersive dielectric response, it is also 



shown that a strong localization occurs for intermediate 
degree of disorder in metallic layers, which was not found 
previously in dielectrics. The findings suggest the possi- 
ble applications over optical devices operating near pho- 
tonic band edges in dispersive periodic-pa-average struc- 
tures in parallel with photonic crystals^. Furthermore, 
the SPS is not rigorously met in the presence of multiple 
scattering effects or residual absorption in our studies. 

The work was supported by National Center for The- 
oretical Sciences (Physics Division) in Hsinchu, Taiwan. 
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